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Abstract 

Using the 9fb~ 1 data sample collected with the CLEO II. V detector at the Cornell Electron 
Storage Ring, we have studied the resonant substructure of the Cabibbo suppressed decay D° — » 
7T~7r + 7r°. We observe significant contributions from the p~ir + , p + ir~, /0°vr°, and non-resonant 
channels, and present preliminary results of the amplitudes, phases, and fit fractions for these 
sub-modes. No evidence for the cr(500) or more massive p resonances was found. We observe no 
CP violation, finding A C p = O.Oljftjj? ± 0.09. 



*Submitted to the International Europhysics Conference on High Energy Physics, July 2003, Aachen 
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I. INTRODUCTION 



A. Motivation 

Three-body decays provide a rich laboratory in which to study the interference between 
intermediate state resonances, and can provide a direct probe of final state interactions. 
When a particle decays into three or more daughters, intermediate resonances dominate the 
decay rate and will cause a non-uniform distribution of events in phase space when analyzed 
using a "Dalitz plot" technique. 0] Since all events of a particular decay mode have the 
same final state, multiple resonances at the same location in phase space will interfere. This 
provides the opportunity to experimentally measure both the amplitudes and phases of the 
intermediate decay channels, which in turn allows us to calculate their fit fractions as well 
as set limits on CP violating asymmetries. 

This paper describes a study of the underlying structure in the decay D° — > 7r - 7r + 7r°, 
and is motivated in part by E791's recent work on D + — > 7r _ 7r + 7r + which found significant 
evidence for a broad neutral scalar resonance, the cr(500).[j| 

Since D° — > 7r~7r + 7r° is a CP eigenstate, we also look for CP violation. With no CP 
violation, D° — > p + 7T~ should have the same amplitudes and phases as D° — > p~ir + . Recent 
theoretical works suggest that CP violation in D° — > 7r~7r + 7r° may be as large as 0.1%. 

a a. 



B. Three Body Decays 

In this analysis we are studying the decay of a spin zero D into three spin zero 7r's, 
hence two degrees of freedom are needed to completely describe the system. 1 A convenient 
choice is to pick any two of the irir invariant mass squared terms, for example m^._ 7r+ and 
m^ +7r o since, when averaged over intermediate spin states, the phase-space for the decay 
is independent of position in (m^.^, m^ +7r0 ) space. A Dalitz plot is simply a scatter plot 
of all event candidates in the {rn 2 K - 1T+1 m^ +7r0 ) plane. Since phase-space is uniform in these 
variables, any structure that shows up in the Dalitz plot is due entirely to the matrix element 
of the decay. Intermediate resonances will show up as bands on the Dalitz plot. Loosely 
stated, the structure within each band tells us about the spin of the resonance, the number 
of events in each tells us about relative amplitudes, and the regions of overlap between the 
bands contains information about relative phases. 

Figure ^ shows a Dalitz plot of all events passing the analysis selection requirements 
in this analysis. It is clear that the distribution of events is not uniform. Indeed, bands 
centered at mass-squared values appropriate for p(770) mesons can clearly be seen along 
both axes. The fact that these bands are not uniform, but rather are populated toward the 
edges of the plot, indicates that they are due to spin-1 resonances. The narrow uniform 
band at m^.- 7r + ~ 0.25 GeV 2 is due to D° — > K° s -k° decays. 



Each daughter has 4 degrees of freedom, for a total of 12. Energy and momentum conservation fixes 4 
of these. Knowledge of the daughter masses fixes 3 more. Since overall spatial orientation is irrelevant, 3 
more are fixed, leaving 2. 
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C. Resonances 

The technique used in this analysis is to fit the (m^- +, m 2 +n0 ) event distribution from 
data with various models containing collections of possible resonant as well as non-resonant 
contributions. Any non-resonant contribution is assumed to be uniform across the Dalitz 
plot, and all resonances are represented by a Breit-Wigner lineshape times appropriate 
factors: 

Mdo^ab)C = F D o(q 2 )F res (q 2 )— ^ (1) 

m res ~Q - im Tes T{q) 

where q 2 = m? AB is the reconstructed invariant mass squared of the (AB) candidate, m res is 
the nominal mass of the resonance, F D o(q 2 ) and F res (q 2 ) are form-factors 0, and T(q) is a 
mass-dependent width. The term in the numerator, As, comes from angular momentum 
considerations and depends on the spin S of the resonance. For spin resonances, Aq = 1, 
and for spin 1, A\ — m 2 AC — m 2 BC + {m 2 D0 — m 2 c ){m 2 B — m 2 A ) /m 2 es . 

It is unknown, a priori, how much of each possible resonance will be required to fit the 
data, thus we weigh each Breit-Wigner term in the matrix element with its own amplitude 
A and phase <fi. The total matrix element is simply a sum of the matrix elements for all 
resonances under consideration in a given model: 
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jvijjv — >-t^ 7r^"7T u ^non res ^ ° ~ 

^p+tt- x *" + ^ x Af230_ p+ir - + (2 
Ap-^+ x e* x A^oo^-^ + 1 j 

^IpO^o X ^f - X A< D°-+p°n° + • • • 

It is clear from a simple examination of the Dalitz plot shown in Figure Q that there are 
at least three resonant contributions to the D° — > 7r~7r + 7T° decay: p°7r°, p + 7i~, and p~7i + . 
Not obvious is the extent to which other resonances may also be contributing, and the level 
at which a non-resonant component is present. As described below, our procedure was to 
start with the three ptr resonances and add other components one at a time to see which 
improved the fit and yielded fit fractions which were not consistent with zero. 

The long lifetime of the K s means that the two-body decay D° — > Kg-K will not interfere 
with any other resonances, thus we treat the Kg as part of the background. 



II. ANALYSIS 

A. General Overview 

We are searching for the following decay chain: 



D* + - i>7T slow 



o 

" slow 





1 — ► 7T 7T ' 7T 

' > 'V 'V 

We require that the D° mesons under study come from decaying D* +, s. This additional 
constraint greatly reduces the background without a large signal loss, and since the charge 
of it^ low is correlated to the charge of the charm quark in the D meson, it also provide a 
method for differentiating between D° — > 7r~7r + 7r° and D° —> 7T + 7r _ 7r° which are otherwise 
indist inguishable . 

We first construct tt° candidates from photon candidates in the electromagnetic calorime- 
ter. Next, we construct D° candidates by combining 7r° with pairs of oppositely charged 
tracks. D* +, s are built by adding 7rf low candidates to the D 0, s. Combinatoric background is 
significantly reduced by demanding that the mass of the D and the D* — D mass difference 
both be within one standard deviation of the nominal value, and that the scaled momentum 
of the D*, x D * = PD*/PD*,max, is greater than 0.7. 



B. Dalitz Fitter 

To fit the data to the matrix element model shown in Equation El we used a MINUIT- 
based unbinned maximum likelihood fitter, minimizing 

T= £ (-21n£), (3) 

events 



Charge conjugation implied throughout. 
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where 



•^"signal ■^'background 



represents the likelihood that a given candidate is either signal or background as a function 
of the fit parameters which determine M.£>o_, n - n+n o, an d 

• F is the fraction of signal events in the sample, about 80% in this analysis, obtained 
by fitting the D° — > 7r~7r + 7T° mass distribution. 

• £^(m^_ 7r+ , m^. +7r o) is the efficiency for an event falling at point (m^ n+ , m^ +7r0 ) in the 
Dalitz plot to be detected by CLEO and to pass all of our analysis cuts. This shape 
was determined by fitting a 2D cubic polynomial to reconstructed signal Monte Carlo 
generated uniformly in phase-space. 

• B{m^- 1T+) m^ +7r o) is the background level at point (ml- n+ , m^. +7r0 ). The background 
shape was studied using data from a D* + — D° mass difference sideband, and was 
parameterized by a cubic polynomial plus additional terms representing real Kg and 
p meson decays. 

• Nsignai = f£(ml- w+ , m£ +7r o) | M D o^ n - n+n o \ 2 dVP is the signal normalization. 

• ^background = JB(ml- n+ , m^ +7r o )dVP is the background normalization. 

III. PRELIMINARY RESULTS 

The preliminary results shown below represent the full CLEO II. V 0, 0] dataset of 9 
fb~ x . Figure 121 shows the three binned Dalitz plot projections for both the data and the 
best fit. 

The only significant contribution to the resonant substructure of this decay was seen from 
the D° — > p + 7r~, D° — > p~7r + , and D° — > p°n° channels. The preliminary fit results are 
presented in Table [I] 



Resonance 


Amplitude 


Phase(°) 


Fit Fraction(%) 


P + 


1. (fixed) 


0. (fixed) 


76.5 ±1.8 ±4.8 


P° 


0.56 ± 0.02 ± 0.07 


10 ± 3 ± 3 


23.9 ±1.8 ±4.6 


P 


0.65 ± 0.03 ± 0.04 


-4±3±4 


32.3 ±2.1 ±2.2 


non res. 


1.03 ±0.17 ±0.31 


77±8± 11 


2.7 ± 0.9 ± 1.7 



TABLE I: Preliminary Fit Results 

Adding other resonances to the fit, including a scalar <r(500), did not result in a signifi- 
cantly improved likelihood and yielded fit fractions consistent with zero. 3 



3 Work to produce upper limits is under way. 
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In addition, there was no evidence of CP asymmetry, either by comparison of the phases 



and amplitudes from separate fits to our D° - 
or in the single calculated figure of merit Acp 



TV TV TV 

-- 0.01±g;g? ± 0.09 



and D° — > tv + tv~ 
where 



TV data samples, 



A 



CP 



\M D o\ 



\M 



DO 



M 



rdVP 



(WP. 



(5) 



We estimate systematic errors by varying parameters of concern and re-running the Dalitz 
fitter to assess the sensitivity of the results to these variations. This was done for our 
parameterizations of efficiency and background, our determination of signal fraction, and 
our event selection criteria. Further systematic studies are ongoing. 



IV. CONCLUSIONS 

We report preliminary results from our Dalitz analysis of D° — > tv~tv + tv°. While we 
find large fit fractions for the three p(770) channels D° — > p + tv~, D° — > p~tv + , and D° — > 
p°7v°, as well as a small but significant non-resonant contribution, we do not find significant 
contributions from any other resonances including the a(500). 

The <t(500) remains a controversial particle. E791 found strong evidence for it in their 
D + — > tv~tv + tv + analysis^, while we find no need for it in D° — > 7r~7r + 7r°. This clearly 
demands further study. 

The lack of a contribution from more massive p mesons is also worthy of comment. In 
CLEO's analysis D° — ► A'~7r + 7r°,|3] the p + (1700) had a small yet statistically significant 
fit fraction. Having observed D° — > A'~p + (1700), one might expect to observe D° — > 
7r~p + (1700) as well, but we do not. It is interesting to note that in the D° — > K~tv + tv° 
analysis, the nominal peak location of the p + (1700) was not on the Dalitz plot so this 
analysis was only sensitive to the tail of this resonance. In D° — > tv~tv + tv°, the peak of the 
p(1700) resonances are located within the Dalitz plot, making D° — > tv~tv + tv° a potentially 
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more suitable laboratory for studying these heavy resonances. Again, this begs for further 
investigation. 

Finally, we saw no evidence for CP violation, either from comparing the amplitudes, 
phases, and fit fractions from the separate D° —> 7r~7r + 7r° and D° — ¥ fits, or in 

the single Acp number. Recent theoretical work permits CP violation on the order of 0.1% 
0, IH in this mode, but we do not have sufficient statistics to confront this prediction. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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